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Abstract 

We exploit the grassmannian nature of the variables involved in the path integral expres- 
sion of the grand canonical partition function for self-interacting fermionic models to show, in 
one-space dimension, a general relation among the terms of it expansion in the high tempera- 
ture limit and a combination of co-factors of a suitable matrix with commuting entries. As an 
application, we apply this framework to calculate the exact coefficients, up to order /3 3 , of the 
expansion of the grand canonical partition function for the Hubbard model in d = (1 + 1) in 
the high temperature limit. The results are valid for any set of parameters that characterize 
the model. 
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1. Introduction 

A quantum system at thermal equilibrium can be completely described provided that one 
knows its grand canonical partition function, which can be expressed as a path integral. For 
bosonic systems, an advantageous feature of the path integral approach is that of employing 
commuting functions instead of non-commuting operators. For fermionic systems, however, 
such an advantage is not obvious to hold, as the integration variables are also non-commuting. 

In Condensed Matter Physics we have self-interacting fermionic models that describe 
correlated electrons. Certainly one of the most studied models is the Hubbard model [1] in a 
variety of space dimensions. 

Due to the non-commutative nature of the fermion fields, it is a common practice to 
bosonize the fermionic model. In the path integral formulation of the Hubbard model, this 
path integral can be re-written in terms of auxiliary bosonic fields using the Hubbard- 
Stratonovich transformation [2]. The Hubbard-Stratonovich transformation allows the in- 
troduction of different decompositions through different auxiliary bosonic fields [3]. Such 
decompositions are equivalent only when the path integral can be calculated exactly. How- 
ever, due to the quartic interaction term in the Hubbard model, only perturbative results are 
attained, and we can no longer relate terms of the perturbative expansion obtained through 
different decompositions [3]. The ambiguity thus arised is handled by a suitable choice of 
decompositions, guided by the physical symmetries of the problem under consideration. 

In the recent literature, we find many papers discussing the high temperature expansion 
of the Hubbard model in space dimensions larger than one [4]. Even though we have exact 
results for the Hubbard model in one space dimension at zero temperature [5], the study of 
the behavior of the unidimensional Hubbard model in the high temperature limit, we only 
find in the work by Takahashi [6] and Shiba-Pincus [7] in the seventies. 

Takahashi [6] derives an integral equation for the thermodynamic potential. However, 
closed expressions for the physical quantities are only obtained in two limits: U — > and 
U — > oo. Shiba-Pincus [7] made a numerical analysis for the unidimensional half-filled- 
band Hubbard model as a function of temperature. They considered two kinds of boundary 
condition: i) finite chain with open ends; ii) finite ring with cyclic boundary conditions. In 
both cases, the model studied has at most six atoms. The results of the finite size model were 
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used to extrapolate to the case of an infinity chain at finite temperature. 

We exploited in a previous paper the grassmannian nature of the fermionic fields to study 
the high temperature behavior of the grand canonical partition function for the anharmonic 
fermionic oscillator [8]. That model has zero (space) dimension, and a natural extension is 
the unidimensional Hubbard model [1]. 

In a recent paper [9], we calculated the exact first-order coefficient in (3 {(3 = -X,) of 
the grand canonical partition funtion for the unidimensional Hubbard model in the high 
temperature limit. Our result was free from any ambiguity, like the one introduced by the 
Hubbard-Stratonovich transformation [3]. Besides, it was valid for any set of parameters 
characterizing the model. 

In this article, the method for evaluating the path integral of self-interacting fermionic 
systems presented in references [8,9] is improved, being thus extended to get the exact higher 
order coefficients of the grand canonical partition function for the unidimensional systems. 

Exemplifying the application of our method, we consider the unidimensional periodic 
Hubbard model. The exact coefficients of order f3 2 and (3 s of the grand canonical partition 
function expansion are calculated in the high temperature limit. 

In section 2 we present the method of evaluating multivariable fermionic integrals, based 
on the grassmannian nature of the fermionic fields. The results derived in this section are 
valid for any unidimensional self-interacting fermionic model. We also write down the kernel 
of the Grassmann function associated to K in the Hubbard model, where K = H — /xN (see eq. 
(2.2)). In section 3 we discuss the symmetries in one-space dimension that allow us to show 
that many multivariable integrals involved in the calculation are equal. Those symmetries 
guarantee that we have a small number of integrals to calculate. In section 4 we calculate 
the exact coefficients of order (3 2 and (3 s in the expansion of the grand canonical partition 
function in the high temperature limit. In section 5 we obtain the Helmholtz free energy of 
the undimensional Hubbard model up to order /3 3 . The result is valid for any number TV of 
space sites. We derive the average energy by site, the average diference of number of spins 
up and down per site and finally the average of the square of magnetization per site. In 
Appendix A, we show that the moments of gaussian Grassmann multivariable integrals are 
equal to co-factors of an appropriate matrix. 
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2. Expansion in the High Temperature Limit and the Grassmann Multivariable 
Integrals 

We have recently shown that the grand canonical partition function of the Hubbard 
model in d = (1 + 1), up to first order in f3, for a finite number iV of space sites is [9], 

„ N M-1 

z iPiv) = / n n n ^{xi^d^xi,^) 

J 1=1 <j=±l i=0 

N M — 1 M— 1 

E E E 4>< T (x l ,T x )(ip a (x l ,T i )-4> a (x l ,T i+1 )) - e K(i>a(xi,Ti),ip a (xi,Ti)) 

ei=lCT=±li=0 e i=0 +0((3 2 ), (2.1) 

where (3 = , /c is the Boltzmann constant, and T is the absolute temperature. The space 
lattice has N sites and M is the number of sites in the temperature lattice, such that eM = (3. 
K is given by 

K = H-^N, (2.2) 

H is the hamiltonian of the system, /j, is the chemical potential and N is the total number of 
particles operator. 

Using a previous result [10], we were ible, for fixed N mnd M, to re-write eq.(2.1) as a 
sum of determinants. We used the simbolic manipulation language MAPLE V.3 to calculate 
the determinants for fixed numerical values of iV and M, and get a recursive expression for 
the grand canonical partition function of the unidimensional Hubbard model up to order (3: 



2N 



l-N0((Eo-p) + j) 



(2.3) 



It is important to point out that the evaluation of determinants of NM x NM matrices 
were in order. Due to memory and speed limitations, such task got restricted to matrices for 
which NM < 30. It was also a matter of luck to obtain the expression (2.3). 

Now, we present a new method to obtain the exact coefficients of the expansion of the 
grand canonical partition function for any self-interacting fermionic unidimensional model in 
the high temperature limit. The fermionic models have iV space sites, where N is any integer 
value. From the results derived we can calculate the thermodynamic limit of iV — > oo. 



The grand canonical partition function in the limit of high temperature is, 



= Tr[l - pK] + ^-Tr[K 2 ] - |-Tr[K 3 ] + • • • , (2.4) 

where K is given by (2.2). Really, (the second line on the r.h.s. of eq.(2.4) is the expansion of 
the exponential of an operator. 

The trace of any fermionic operator, written in terms of creation (aj) and destruction 
(a_j) operators, can be mapped in terms of the generators of the Grassmann algebra {fji,rjj}, 
provided we make the identification: 



aj -> r]i and a.j -> — . (2.5) 

Let us consider a Grassmann algebra of dimension 2 2Af , whose generators are: {fji, ■ ■ ■ , 
W; Vii " ' ' > w}- The generators rji, fjj satisfy anti-commutation relations: 



{Vi, Vj} = 0, {Vi, Vj} = and {rji, rjj} = 0, 

where i, j = 1, 2, • • ■ , Af. 

The product of two fermionic operators A and B, written in terms of the Grassmann 
generators, is 

M 

(AB)(f), v )= / Yldri'jdfj'jAfav')^'^) e ^ , (2.6) 
i=i 

where A(fj,rj') and B(fj',r}) are the kernel of the fermionic operators A and B, respectively. 
For fermionic operators in normal order, 



A = E (*\) ni ■ ■ ■ (4f) n " M m "---M mi , (2.7a) 

and 

B = E B t,-,n N ; mi ,..;m N (4)™ ' ' " ' ' ' , (2-76) 



rn^ , ■ ■ ■ , ?Ti^y- =0 



l 



«l.---.«A/'= 

771 ]_)■■■ ,rnj\f=§ 
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with A® 5 ... )rw;mi5 ... )mjv . and B$ u ... jn „- mi ,-,mM being commuting constants. The kernel of 
these operators are given by, 



A(fj,v') = e^ A®(fj,v') (2.7c) 

and 

E v'lVl 

B{f)',ri) = e'=i B®(fj / ,rj) (2.7d) 

where, in a naive way, we can say that we get A®(fj, rj') from eq.(2.7a) by replacing aj — > f/j 
and aj — > r/j, and similarly for B®{f)',rj). 

All the operators considered in this paper are in normal order. 

The trace of any fermionic operator O is [11] 

AT 

YldriidfjiOfari) e'=* , (2.8) 
i=i 

where we use the shorthand notation: fj = {fj\, • • • , fj_\f} and rj = {771 , • • • , ?7at}, and 0(fj, rj) is 
the kernel of the fermionic operator O. For the case where the operator O is a product of n 
normal ordered fermionic operators Q, we use the relations (2.6-8) and algebraic manipula- 
tions to write the trace of such product of operators as 



r^L n ^} EE Vi(t„)[vi(t„)-Vi(t„+i)] 

TMQ n ] = / II II dviMdfjii^) e«=i -=° x 

1=1 n=0 

x O®(f)(r ),r 1 (T ))O®(f}(Ti),v(ri)) x ... x 0®(f)(r n _ 1 ), V (T n _ 1 )), (2.9) 

where we use the convention: fj(r u ) = {vi(r u ), • • • , fjtf(T u )} and rj(r u ) = {?7i(tv), • • • , w( r ^)} 
for v = 0, 1, • • • , n — 1. The Grassmann variables in eq. (2.9) satisfy the boundary condition: 

Vl( T n) = -Vi( T o) and t]i(t u ) = 0, for 1/ > n, (2.9a) 

and 1 = 1,2, ■■■,M. 
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Relation (2.9) is used to write the terms of the expansion of the grand canonical partition 
function in the high temperature limit as multivariable Grassmann integrals. 

For a one-dimensional fermionic model, the index I of the generators of the Grassmann 
algebra: {t7z(t m ), f)i(r M )}, conveys space and spin degrees of freedom. Explicity separating 
the space index from the spin index in the generators of the non-commuting algebra, the 
coefficients of the expansion Z{(3,[i) in eq.(2.4) become, 

JV n-1 

/N n-1 ^ ^ ^ fj tT (x l ,T u )[r) tT (x l ,T v )-ri< T (xi,T v+1 )] 

nn n ^^' r -)^(^' r -) e!=iCT=±ii ' =o x 

1 = 1 <T = ±1 V=0 

x JC®(f) a (x,To),r] <T (x,To))}C®(fi (T (x,T 1 ) 1 r] <T (x,T 1 )) x ••• x JC®(fj a (x, r n _i), r] a (x, r n _i)), 

(2.10) 

where iV is the number of space site, a =|= +l,cr =|= —1, and the boundary conditions 
(2.9a) become, 

Va(xi,r n ) = -r] a (xi,r ) and r] a (x h t v ) = 0, for v > n, (2.11) 

/ = l,2,---,iV. 

To make use of previous results, where we relate the multivariable Grassmann integrals 
to a sum of determinants [10], it is interesting to map the generators f] a {xi, r„) and fj a (xi, r„) 
into single-indexed anti-commuting variables. The particular mapping is a matter of choice; 
however, calculations are greatly simplified if we choose: 



rft(xi,T v ) = ri vN+l (2.12a) 

and 

rn(xi,T v ) = ri( n+v ) N+l , (2.126) 

where I = 1, 2, • • • , N. N is the number of space sites, and, v = 0, 1, • • • , n — 1. 
The mappings (2.12 a-b) can be summarized as: 



Va(xi,T„) = V[ ( 1 -c [ln+u]N+l 
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(2.12c) 



Using the redefined generators (2.12c), the sum in the exponential on the r.h.s. of 
eq.(2.10) is written as 



n-1 



Yl ^2vA x ^ T v)[V<ri. x h T v)-vA x ^ T v+l)] 
:1 v=0 

2nN 

Y Vi A u Vj, 



1 = 1 cr = ±l v = 
2nN 



(2.13) 



i,j=i 



where Ajj are the entries of the block-matrix A, 



A = 



(2.13a) 



q A u 

whose entries are matrices of dimension nN x nN. The indices /, J are such that I, J = 
l,2,---,nN. 

The matrices A^ and A^ are identical and have a block-structure. Taking into account 
the anti-periodic condition in temperature (2.11) in (2.13), the matrices A <Ta , a =\, \, are 



A^ = A^ 



\ljVxJV QvxJV QvxN 

The non-null elements of A acr , a =j and a =|, are: 



QvxTV \ 
QvxTV 

liVxAf / 



(2.14a) 



a n = 1, / = 1,2, • • -,nN 

AYj = I a IJ+N = -1, /= l,2,..-,(n-2)JV 

, «(n-l)iV+/,7 



(2.146) 



1, 7=l,2,--.,iV 



The matrices A aa , a =t,|, have dimension nN x nN. The symbols 1tv X 7V and GhxN 
stand for the identity and zero matrices of dimension N x N, respectively. 

With the newly indexed generators, the expression of Tr[K n ] (eq.(2.10)) becomes, 



2nN 

/2nN ^ fjj Au rjj 

J I drudfjj e'' J=1 x 
1=1 

x K,®(fj, 77; v = 0) /C®(ij, 77; 1/ = 1) • • • /C®(rj, n;u = n-l). (2.15) 

Note that expression (2.15) is the exact coefficient in order f3 n of the expansion in the 
high-temperature limit of the grand canonical partition function for any unidimensional self- 
interacting fermionic model. The specific model to be studied is represented by the function 
/C®. The matrix A is the same for all unidimensional fermionic models. 

The Grassmann functions K.® are polynomials in the generators of the algebra. Therefore, 
the r.h.s. of eq.(2.15) are moments of the multivariable Grassmann gaussian integrals. We 
show in Appendix A that these integrals can be written as co-factors of the matrix A. 

Once the sub-matrices A^ and A^ are null, the multivariable integral (2.15) is equal 
to the product of the contributions coming from the sectors: aa =f| and aa =||. 

Even restricting ourselves to the calculation of the multivariable integrals of a fixed 
sector aa, we have to obtain the determinant of non-diagonal matrices of dimension nN x 
nN. The evaluation of such determinants by means of some symbolic manipulation language 
obviously depends on hardware and software resources. For a fixed n, for instance, we have 
an upper practical limit for iV for doing the calculation of the determinants. One possibility 
for evaluating (2.15) is to fix different values for N and try to extrapolate the results for an 
arbitrary value of N. If we are lucky, we may recognize some recursion expression for (2.15) 
for all N. 

The integrals in (2.15), for each sector aa, have the form: 

nN 

„ nN ViAijijj 

M(L,K) = / Hdf]^ f) hVkl ---f)i m Vk m e-= 1 , (2.16) 

J i=i 

with L = {Zi, • • • , l m } and K — {ki, ■ ■ ■ , k m }. The products fjrj are ordered in such a way 
that h < fa < • • • < l m and k\ < k 2 < ■ ■ ■ < k m . From eq.(A.lO), the result of this type of 
integrals is equal to: 

M(L, K) = (-i)( l i+h+-+i m )+(ki+k 2 +-+k m ) A ( L ^ R ^ ( 2 _ 17 ) 
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where A(L, K) is the determinant of the matrix obtained from matrix A by deleting the lines 
0i> • • • j lm} and the columns: {ki, • • • , k m }. 

Our approach to calculate the integral (2.16) is, for fixed n and arbritary N, to explore 
the block-structure of matrices A aa , a =\ and a =|, diagonalizing it through a similarity 
transformation 



p- 1 A fJfJ P = D, (2.18a) 

where the matrix D is, 



D = 



(2.186) 



Xi, i = 1, 2, • • • , n, are the eigenvalues of matrices A "' 7 , a =f, |, and calculate the co-factors 
of the matrix D. 

The j t/l -column of matrix P is the eigenvector of A aa associated to the eigenvalue Xj. 
Each eigenvalue of matrix A aa has degeneracy N. The matrices are not hermitian, thus some 
eigenvalues are complex. In the following, we will be working on the ao =|T sector; however, 
the results for the a a =H sector are analogous since A^ = A-'--'-. 

We make a change of anti-commuting variables, 

^^P" 1 ^ and fj' = f]~P, (2.19) 

where rf = {r)[, • • • , r)' nN } and fj' = {fj[, ■ • • , fj' nN }- The jacobian of the transformation (2.19) 
is equal to one. 

In eqs.(2.19), we wrote the transformations of variables in a very simplified way; however, 
due to the fact that A^ is a block-matrix, the matrix P also has a block structure. This fact 
implies that transformations (2.19) do not mix the index associated to the space site. 

In a schemmtic way, the integrals M(L,K) (eq.(2.16)) become: 



„ nN 

M(L,K)= / Hd Vi dfji (r ? P- 1 ); 1 (Pr / ) fcl ---(r ? p- 1 ) ;m (P7 ? ) fcK 

J i=i 



nN 

e^ =1 , (2.20) 

where Dij are the entries of matrix D. 

From eq.(2.16), we have that M(L,K) is a multivariable Grassmann integral, which in 
turn correspond to co-factors of the diagonalized matrix D (eq.(2.17)). It is very simple to 

2 2 

calculate these co-factors. Besides, taking into account that r]' = and fj' = decreases 
the number of integrals to be calculated in (2.20). 

This is a general approach, and it can be applied to any self-interacting unidimensional 
fermionic model. Thus we have reduced the calculation of fermionic path integral to the 
evaluation of co-factors of a diagonalized matrix. 



2.1. Application to Unidimensional Periodic Hubbard Model 

We will apply in section 4 the general approach presented in section 2 to the periodic 
Hubbard model [1] in d = (1 + 1) with iV space sites in the presence of an external magnetic 
field in the z direction. However, we need to know the normal ordered Grassmann function 
K.® for the specific model. 

The hamiltonian that describes the Hubbard model in one space dimension is [1]: 

N N N 

H = 5^ ^> a :/> + ^XXW a *W + Xb Y1 Yl aa t a - ( 2 - L1 ) 

i = l (T=-l,l i=l i = l (T=-l,l 

J = l 

where aj^ is the creation operator of an electron in site i with spin [12] a and a icr is the 
destruction operator of an electron in site i with spin a. All diagonal elements of tij are 
equal, tu = Eo, where Eo is the kinetic energy. The only non-null off-diagonal terms are 
Uj-i = ti : i + i = t, where i = 1, 2, . . . , N, and they contribute to the hopping term. U is the 
strenght of the interaction between the electrons in the same site but with different spins. We 
have defined = -i^B, where g is the Lande's factor and hb is the Bohr's magneton. 

The periodic boundary condition in space is implemented by imposing that a 0(T = &Na 
and a7v + i iCr = a^. Therefore, the hopping terms tioaj^aocr and tjv,iv+ia] Vcr aiv+i,CT become 
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tiAraj^ajvcr and tjv.iajvo- 21 !^ respectively. We point out that the hamiltonian (2.1.1) is already 
in normal order. 

From eq.(2.2) we have 

K = H - /iN. 

The kernel of the operator K of a one-dimensional Hubbard model in a lattice with iV 
space sites, written in terms of the generators r} a (xi,T„) and fj a (xi,r u ) is 



N 



K®(fj a (xi,T v ),ri a (xi,T v )) = ^2^2(E + a\ B - n)fj a (xi,T v )ri a (xi,T v )+ 

1 = 1 cr±l 

JV 

+ t[w(xhT„)r] a (xi +1 ,T„) + fj a (xi,T v )ri a (xi-i,T v )] + 

1 = 1 <T±1 

JV 

+ ^2 u V]ixi,T u )rj ] ixi,T 1 ,)f] l (xi,T iy )r] l (xi,T iy ), (2.1.2) 
1=1 

with the space periodic boundary condition imposed by: 

tfj a (x N , r„)r] a (x N+1 ,T u ) = tf} a (x N , T v )r} a {x-i, t„) (2.1.2a) 

and 

tfj a (x 1 ,T„)ri a (xo, t v ) = tfj a (xi, T„)r) a (x N , t v ), (2.1.26) 
and the anti-periodic boundary condition in r: 

V(r(xhT n ) = -7/ ff (xi,r ), (2.1.2c) 

for I = 1,2, and a =|,|. 

Using the mapping (2.12c), the kernel of the operator K becomes: 

N 

1 = 1 a±l 
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N 

Z=l <T±1 

TV 

+ ^ Ufj(n+v)N+l V(n+v)N+l VvN+l VuN+l- (2.1.3) 



3. Useful Symmetries of the Mult i variable Grassmann Integrals 

In eq.(2.15), the number of terms to be calculated increases rapidly with N. To make 
this approach feasible, we have to explore symmetries, besides the trivial space translation, 
of the exponential in eq.(2.15). In section 2, we pointed out that in expression (2.15) the 
contributions from the sectors aa =]] and a a =H are decoupled. In this section we discuss 
only one of the two sectors, e.g. oo =\\. Analogous results are valid for the sector aa =||. 
In the expression of Tr[K n ] (eq.(2.15)), we deal with integrals of the following type: 



/nN 
drudfji fj^N+hVviN+k! Vv2N+i 2 Vv2N+k 2 x • • • x 
7=1 

nN 

E A\\ nj 

x V» m N+i m V» m N+k m x e'" 7=1 , (3.1) 

where m < n, and, v\ < u 2 < • • • < v m . 

We want to study how the integral I(m) transforms under a change of variables in the 
temperature parameter r. Consider the change of variables: 



vIn+i = V{*-i)N+i, i/ = l,2,---,n-l, (3.2a) 

and 

V'l = -V(n-i)N+h (3-26) 



V'uN+l = V(u-l)N+h 1/ = 1,2, 1, 



(3.2c) 



and 
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v'i = -V(n-i)N+i- (3.2d) 

In this relabelling of generators, we are shifting the temperature index by one unit, except 
for r) n -i and r\ n -\i identified to — rjo and — fjo, respectively. The jacobian of the simultaneous 
transformations (3.2) is equal to one. 

Before discussing how the integral I{m) transforms under (3.2), let us show that the sum 
in the exponentiml of eq.(3.1) i invarimnt under such transformations. 

From eq.(2.14b), we have that the non-null terms of A^T are: 

a H = 1, / = 1, 2, • • • , nN 

ai,i+ N = -l, 1 = 1,2, •••,(n -2)N (3.3) 
, a(n-i)N+i,i = 1, J =1)2, ■■■,N 



Ada 



Therefore, the sum within the exponential in eq.(3.1) can be written as, 



nN n-2 N 



1,3=1 u=0 1=1 

N 



w A l j vj = Yl Yl vvn+i A V N+UvN+i vvn+i+ 

v=0 1=1 
N 

+ Y f )(n-l)N+l ^(I_i )A r +Z)(n _i )A r +Z V(n-l)N+l + 
1 = 1 

n-3 N 

+ YY^ N + l A In+1,(u+1)N+1 V(u+l)N+l + 
u=0 1=1 

N 

+ Y f j(»-2)N+l A ll- 2 )N+l,(v-l)N+l r l(v-l)N+l + 
1=1 
N 

+ J2v(n-i)N+i A\l_ 1)N+ljl m. (3.4) 



1=1 

Under the change of variables (3.2), the terms on the r.h.s. of (3.4) transforms as: 

n-2 N 

1 )J2J2^ N + l A lN+l,uN+lVuN+l = 
v=0 1=1 

n-2 N 

VvN+l A ( u -i)N+l,(is-l)N+l VuN+l 



v=0 1=1 
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n-1 N 



^2^2 V'uN+l ^vN+l,uN+l WvN+l- 



v=l 1=1 



N 



r7jv( n _i) +z A N{n _ 1)+ ^ N{n _ 1)+l r) N ( n _ 1)+ i - 

1=1 

N 

= E * 4! 



-3 AT 



3 ) XI Yl VvN+l A Vn+1,{u+1)N+1 7] {l/+1)N+l = 
u=0 1=1 

n-1 N 

= EE 

VvN+l A {v-l)N+l,vN+l V( v +l)N+l 

v=l 1=1 
n-1 N 

= ^2^2 VuN+l A uN+l,(u+l)N+l V[u+1)N+1- 



i=l 1=1 



N 



El T 
V(n-2)N+i A (n _ 2)N+l ^ n _ 1)N+l r] {n _ 1)N+ i = 

1=1 

N 

= _ 5Z A \n-l)N+l,l Vn+1 
1=1 

N 

= X Vi a JJn+i Vn+i ■ 



i=i 



N 



5 ) V{n-i)N+i ^ n - 1)N +i,i m = 

1=1 

N 

= ~~ X/ V(n-l)N+l A (n-1)N+l,(n-l)N+l Vl 
1=1 

N 

= X V(n-l)N+l A (n-l)N+l,l Vl- 



1=1 
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To obtain the results (3.4a-e), we have used eq.(3.3). 
Summing the r.h.s. of (3.4a-e), we finally conclude that 

nN nN 

E * A u VJ=J2v'i A\l V'j. (3-5) 
i,j=i i,j=i 

We point out that the invariance of the sum (3.5) is true for any value of n and N. 
Moreover, since = A-W, a similar invariance under (3.2) can be written for A-^. 

In (3.2) we translated the temperature index by one unit. However, the equality (3.5) is 
still valid if we translate the temperature index by any integer uq, either positive or negative, 
i.e., 

VvN+i — > r]^ +Uo ) N+ i (3.6a) 

and 

r} 9N+ i -> -rji , where V + v = n, (3.66) 

and perform an analogous transformation for fj. Transformation (3.6) is equivalent to applying 
uq distinct and consecutive transformations (3.2). 

Under the transformation (3.6), the integral I(m) (eq.(3.1)) becomes, 

/nN 
Yl drj'jdfj'j fj[ Ui+Vo)N+h rj[ Ui+Uo)N+ki x • • • x 
i=i 

nN 

E Vi A Vj v'j 

X V(v m+ v )N+l m V(v m+ v )N+k m X e '"' =1 > ( 3 - 7 ) 

where u is a fixed integer, positive or negative; for v such that v + z/ = 0, we assume that 



V(v+v )n+i -»• -m and V(p+u )n+i -»• 

and for v + z/ > n 

V(v+u )N+l ^> V(v+v -n)N+l and 7]( ;y+;yo )Ar + ; — > rf( v+Vo - n )N+l- 
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An analogous result is valid for the sector a a =\\. 

As the contributions to the integrals in eq.(2.15) of the sectors aa =|T an d ccr =|| 
f actor ize, we can choose differents values for i/q f° r each spin sector. 

Another symmetry of the exponential in eq.(3.1) results from the arbitrariness of either 
clockwise or counterclockwise labelling of space sites in the periodic unidimensional chain. The 
result of the integrals is choice-independent. Let us label the space sites counterclockwise. For 
a fixed the Grassmann generators associated to the sites are: {ry^jv+i, 77^+2, • • • ,VvN+n] 
fjvN+i, VvN+2, ■ ■ ■ , VvN+n}- If, on the other hand, we label the unidimensional chain clock- 
wise, in this case the Grassmann generators will be: Wvn+iiVIn+2i ' ' ' iVIn+n'iVIn+h 
vIn+2i ' " 1 vIn+n}- 

The correspondence among generators in the two labellings may be chosen to be 

VvN+i = v'uN+i and VvN+i = f? vN+1 , (3.8a) 
for v = 0, 1, • • • , n — 1, and 

VuN+l = v'uN+(N+2-l) and VuN+l = v'uN+(N+2-l)- (3.8b) 

The jacobian associated to (3.8) is equal to one. 

The sum within the exponential in the integrand of eq.(3.1), namely, 

nN n—1 N 

E VI A Vj ^ = EE VvN+l AIn+1,uN+1 VvN+l + 
I,J=1 v=0 1=1 

n-2 N 

TlvN+l ^-i/N+l,(i/+l)N+l V(v+l)N+l~r 

v=Q l-l 
N 

+ E V(n-i)N+i A\l_ 1)N+ljl vi, (3.9) 
1=1 

is invariant under the substitution (3.8). In order to show this, we use the relation among 
the non-zero elements of A^T, given by eq.(3.3), and manipulate the expressions in a similar 
way as we have done before. In summary, we get that 

nN nN 

E * A u V = E # A ^ Vj, (3-10) 

I,J=1 I,J=1 
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where the relation between {77, 77} and {//, rj'} is given by (3.8). 

The integral I(m) (see eq.(3.1)) subjected to the change of variables (3.8) becomes, 



/nN 
W drjjdfjj ^ lJV+ (j V+2 _/ 1 )^ 1 iv+(Ar+2-fei) x ' ' ' x 
7=1 

nN 

E Vi A Vj Vj 

X V l y m N+(N+2-l m )V l y m N+(N+2-k m ) X e '' J=1 ■ (3.11) 

The chiral space symmetry, represented by the transformations (3.8), is valid for any 
positive integers integers n and N. 
From eq.(2.15), we have that 



/2nN Vi Au r)j 

1=1 

x K®% 77; ^ = 0) /C®(t7, 77; 1/ = 1) ••• /C®(?7,77;z/ = n-l), (3.12) 

where the matrix A is given by (2.13a). For the one-dimensional Hubbard model with N 
space sites, the kernel /C®(r/, 77; v) is given by eq.(2.1.3). 

In order to simplify the algebraic manipulation for calculating eq.(3.12), we define: 



N 

£(fj, 77; v- a) = J2 V { (i^l n+1/]N+l V { i±^l n+v]N+l 5 (3.13a) 
1=1 

N 

T-{n,v;v;<j) = Y,v [ (i^n+v]N + i r i [ (i^i n+u]N+l+1 ; (3-136) 
1=1 



and 



N 



T+(fj, V, ^^) = Y1 ^Ll^ln+^N+t V^^N+L! ■ (3-13c) 



2 

1 = 1 



We call 



S(rj,T];iy) = £ '(a) £(77,77; 1/; a), 



(3.14o) 



r=±l 
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T- (fj, 77; 1/) = t T ~ (V, W v\ °) , (3-146) 

a±l 

T+(fj, m u) = J2t T + (V, V, v\ <r), (3.14c) 



<7±1 

and 

JV 

U(fj,r];u) = ^f] {ri+u)N+ i r] {n+u)N+ i fj vN+ i rj uN+ i. (3.14d) 

where -E(cr) = E — crAs — /x. The term ^(77, 77; v) represents the kinetic energy, T~(fj, 77; v) 
and T + (rj,r]; v) are the hopping terms and U(fj,r]; v) is the fermionic interaction term. 

Using the definitions (3.13) and (3.14), the Grassmann function K,®(fj,rj; v) is written as, 

K®% m y) = £% m v) + r~% v; v) + t+% m u) + u% m (3.15) 

Since the terms in £(77, 77; v), / T~(fj, 77; u), T + (fj, 77; v) and U(fj, 77; v) are products of 77 and 
77 at the temperature index v these expressions still have the same form under the change of 
variables (3.6), except that they are defined at v + vq\ that is, 



£(77,77;*/) -> £(77,77;!/ + 1/ ) (3.16a) 

T" (77, 77; 1/) -> T"(t7,?7;z/ + z/ ) (3.166) 

T+(f),r);v) -> T+ (77,77;^+^) (3.16c) 

U(fj,T];u) ->• U(fj,r]; v + u ). (3.16d) 

Let us see the behavior of the expressions £(77, 77; 1/), T~(fj, 77; z/), T + (fj, 77; z/) and 
^(77, 77; 1/) under the chiral transformations (3.8). We begin by considering £(77,7751/). From 
(3.13a) and (3.14a), we have that 

£(77,77;!/)= ]T S( ( 7)77 [ ( i -^) n+i/]JV+1 7J [i i^ n+l/]jv+1 + 
cr = ±l 

JV 

<7 = ±1 Z = 2 

The transformation (3.8) for a sector <rcr is, 
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\^l n+v]N+1 - v[(i^) n+iy]N+1 and V { (i^l n+u]N+1 - V[(i^) n+1/]N+1 , 

(3.18a) 

for / = 1 and v = 0, 1, • • • , n — 1, and 

V[ii^l n+v]N+l = v' [ (i^i n+v]N+N+2 _ l and ^o^)^^ = ^(i^) n+I/]Ar+JV+2 _ / , 

(3.186) 

for / 7^ 1 and z/ = 0, 1, • • • , n — 1. 

Substituting (3.18) in eq.(3.17), and defining the index k = N + 2 — / on the second term 
in the r.h.s. of eq.(3.17), we obtain 

£(fj,rr,v)= E ^)\^ n+v]N+ ii^i n+v]N+1 + 

<7 = ±1 

N 

+ Yl J2 E ( (T ) f i{(l^ln+ I ,]N+k T] [^ln+u]N+k ■ 
a=±lk=2 

= £(v',v';v) (3.19) 

£(77,77; v) is invariant under the chiral transformation (3.18). 
Proceeding in a similar way, we also get that 

U(fj, V ;u) = U(fj',v'^), (3.20) 

where the relation between the generators {77,77} and {r)',fj'} is given by (3.18). 
Let us now take under consideration one of the hopping terms, namely, 



T (77,77;*/)= t V [ il_?l n+ „ ]N+l V [ il_?l n+ „ ]N+2 + 

a=±l 

N-l 

+ Yl J2 tf j[(^ln+v]N+iy-±^ln+v]N+l+l + 



a=±l 1=2 



+ t \^ln+ V ]N+N V { il^l n+u]N+N+1 , (3.21) 



r=±l 
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where we have already taken into account the periodic space boundary condition. 

Applying the change of variables (3.18) in eq.(3.21) and defining the index k = N + 2 — /, 
we have 

N 

T-(fj, V ;v) = J2 y^l^M+k 77[(i^) n+I/]JV+fc _i 

fc=l cr = ±l 

= T+(fj, m v) (3.22) 

In a similar way, it may be shown that, under (3.18), T + (f/,?7;i/) transforms into 
T-(fj,rj; v). 

We point out that the results of this section are valid for any positive integer values of n 
and N. 

4. Calculation of the Coefficients of the Grand Canonical Partition Function for 
the Hubbard Model in the High Temperature Limit 

Up to now, all that we have said about Tr[K n ] in sections 2 and 3, is valid for any value 
of n and N. 

The calculation of Tr[K n ], given by eq.(2.15), for arbitrary n can formally be done using 
the result (2.20). However, due to the fact that the number of terms that contribute to Tr[K n ] 
increases rapidly with n, we are able to calculate the exact coefficients up to order (3 3 of the 
grand canonical partition function of unidimensional Hubbard model in the high temperature 
limit. 

When we say that our result is exact up to order /3 3 , we mean that it is valid for any 
value of the constant parameters Eo,t,U and fj,, that characterize the model, and, for any 
value of the constant external magnetic field B. 

4.1. The Exact (3 2 Coefficient in d=(l+l) Hubbard Model 

From eq.(2.15), the expression of Tr[K 2 ] in terms of the Grassmann generators is 
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p 4-/V rn Au r/j 

Tr[K 2 ] = / j'l drudfj! e'- J=1 x K®% 77; v = 0) /C®(?), 7/; 1/ = 1) , (4.1.1) 

^ j=i 

where the matrix A is given by (2.13a), and for the unidimensional periodic Hubbard model, 
K,® is given by eq.(2.1.3) or eq.(3.15). 

For n = 2, the matrices and A^ (eq.(2.14a)) are: 

(^-NxN —^-NxN \ 
(4.1.2) 
InxN ^NxN / 

where 1at x at is the identity matrix of dimension N x N. For the case n = 2, the eigenvalues 
of the matrices A^ and A^ are: 

\ l = X* = l-e^. (4.1.2a) 

The eigenvalues Ai and A2 are iV-fold degenerated. 

^From eq.(2.18a), P is the matrix that, through a similarity transformation, diagonalizes 
the matrices A^ and A.**. For n = 2, the matrix P is, 



and its inverse is: 



IN x N l NxN 

HnxN —i^-NxN 

Itvxtv — Hnxn 



(4.1.26) 



P" 1 = - I I • (4.1.2c) 



1 

InxN iljVxA 

Since P and P _1 are block-matrices, their entries can be written as P U N+i,v'N+k 
>-i 



SikPuv' and P vN+ i u 'N + k = SlkQuv'- Vw' and q uu i are the non-null elements and for n = 2 
they are equal to 



( '• _M and V = d j / ) (1.1.2f/) 



Therefore, 
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(1 - e*2 ) InxN %x7V 



, for a=t,|. (4.1.2e) 



QvxTV (1 - e 2 ) l^x^ 

For n = 2, the change of variables (2.19): r] = Pr/ and fj = rfP -1 , becomes, 



VvN+i = *W' vl'N+i (4.1.3a) 

i>'=0 

and 

l 

fj v N+l = Yl rfv'N+l Qv'v , (4.1.36) 

where v = 0, 1 and Z = 1, 2, • • • , N. 

We should note that due to the block structure of the matrices P and P -1 , the space 
index / in the generators r\ and fj are not mixed up by the transformation (4.1.3). 

From eq.(3.15), we write function K,®(fj, rj; v) as 

K®(fj,T);v) =£(f},ri;v) + T~ (fj, 77; v) + T+ (fj, 77; u) +U(fj,r];u). (4.1.4) 
In order to simplify the notation, we define 



7b I V 

/2nJV ^ ?7j Ajj 77,7 



1=1 

x Ox(fj,7];vx)---O m (f],r];v m ) (4.1.5a) 



and 



/(3-a)nN ^ r?(i-<r)„;v+/ 77 ( i_ (T )„ J v + j 

I I dmdfji e 1 ^- 1 x 

I=(l-a)nN+l 

x O x (fj,r];v 1 ) ■ ■ - O m (fj,r];v m ) (4.1.56) 
where Oj(fj,r\; uj) are Grassmann functions. 
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When we substitute eq.(4.1.4) for v = and v = 1 in eq.(4.1.1), we get sixteen terms. 
However, using the results (3.19), (3.20) and (3.22), and the fact that the terms: 

< £(0)T-(1) >,< 5(0)T+(1) >,< T-(0)T-(1) >,< T-(0)W(1) >,< T+(0)T+(1) >, and 

< T + (0)W(1) > are null, we have that Tr[K 2 ] reduces to the sum of four distinct terms: 

Tr[K 2 ] =< 5(0)5(1) > +2 < T-(0)T + (1) > +2 < 5(0)W(1) > + < W(0)W(1) > . (4.1.6) 

Before calculating these terms, lets us consider one of the null integrals and show the reasoning 
for its vanishing. For example, the term 



< 5(0)T-(1) >= 2t(E - n)[< 5(T, 0)T"(T, 1) > + < 5(T, 0)T"(|, 1) >]. (4.1.7) 

In expression (4.1.7) we identify the following types of terms: 
i) Terms coming from < 5(t,0)T~(|, 1) >: 



p 2N fj! A TJ VJ 

Zi = / Yi dr 1idfji e 7 " 7=1 ViiViiVN+hm+h+i x 

J i=i 

4N 

4N Vi A u V i 



'fj drudfj! e I ' J = 2N+1 , (4.1.7a) 



x 

I=2N+1 



where l\ and li are fixed and represent the space site index. 

Making the change of variables (4.1.3) to diagonalize A' 70 ", a =f, |, we get 



l 

Ji = det(A u ) QviOPouiVviiPiu. 

"1,^=0 



/ X 

2 



2N 

2N £ fj'j Du v'j 

X ' 



/ / j 'ii '/j 
II d Vidv'i e'' J=1 V^N+hV'^N+hVlzN+hV^N+h+i- (4.1.76) 

7=1 

The presence of 77^ and r\j in the integrand of eq.(4.17b) implies that the integral is equal 
to the determinant of a matrix obtained from matrix D by cutting the line % and the column 
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j. As long as D is diagonal, the only non-vanishing integrals are those for which i = j. In 
other words, if the i-th line of D is cut, so must be its i-th column. 

In eq.(4.1.7b) the space index (I2 + 1) never equals any other space index. Therefore, the 
minor obtained from matrix D by cutting the lines: {v\ + 1±) and {yi + 12), and the columns: 
(u[ + li) and (u' 2 + l 2 + 1), will always have, at least, one line or row of zeros. Therefore, its 
determinant vanishes. 

Then, 

X x = =^< £(T,0)T-(T,1) >= 0. (4.1.8a) 
ii) Terms coming from < £(T,0)T-(j, 1) >: 

Using a similar reasoning as for the terms in < £(f, 0)T~(f , 1) > , we get that 

<T"(|,1) > 1= 0. (4.1.86) 
Therefore, using the results (4.1.8a-b), we finally obtain, 

<5(0)T-(1)>=0. (4.1.9) 
Proceeding in an analogous way, it is straightforward to show that the terms 

< £(o)r+(i) >,< r-(o)r-(i) >,< r-(o)w(i) >,< r+(o)r+(i) > and < r+(o)w(i) > 

also vanish. 

The results (4.1.8a-b) are a direct consequence to the fact that the change of variables 
(4.1.3) does not mix up space indices. 

Taking into account the previous results and calculating the coefficients of the terms of 
eq.(4.1.6), Tr[K 2 ] may be written as: 



Tr[K 2 ] = 2((£? - y) 2 + X 2 B ) < £(T, 0)5(T, 1) > +2({E - yf - X%) < £{[, 0) >< £(t 0) > + 
+ At 2 < T~ (T, 0)T+(t, 1) > +4(^o - fj) < £(T, 0)W(1) > + < W(0)W(1) > . 

(4.1.10) 
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The terms on the r.h.s. of eq.(4.1.10) are sums of multivariable Grassmann integrals that 
reduce to the following types: 



f 2N Vi AH VJ 

Wi(0 = J J] dr li df ii e 7 " 7=1 Vi Vl, (4.1.11a) 



7=1 



r 2N E Vi AH VJ 

H 2 (l,k)= / Yl dVidfji e 7 " 7=1 ViVi VN+k VN+k, (4.1.116) 

J T — 1 



7=1 



and 



r 2N E vi AH vj 

H 3 (l,k)= / Yl drfidfji e 7 - J=1 Vi Vi+i VN+k VN+k-i, (4.1.11c) 

J 7=1 

in the sector a a =tt? an d equivalent integrals in the sector J, J,. 

Let us calculate in detail these integrals. Those integrals that contribute to order (3 3 
of the expansion of the grand canonical partition function in the high temperature limit are 
handled in a similar way. In the next section, we will restrict ourselves to give the results of 
such integrals. 

Under the change of variables (4.1.3), the integral Tli(l) becomes: 

2JV 

1 , 27V Vi Du v'j 

w i(0 = Yl VvoPov> / JJ drfidfji e 7 - J=1 ff vN+l vl'N+i- (4.1.12) 

v,v'=0 J 7=1 

We only get non-null results for eq.(4.1.12) when the «-th line and the i-th column of the 
matrix D are cut simultaneously. This only happens when v = v' . 
Therefore, 



2JV 

r 2N E Vi D 'J Vj 

= / dq'jdffj e J - J=1 (gooPoo ViVi + QioPoi Vn+iVn+i)- (4.1.13) 

J 7=1 

The first term in the integrand corresponds to cut the /-th line and l-ih column of matrix D. 
Since it is a diagonal block-matrix, the result of the integral is independent of / and it is equal 
to: Aj v_1 A^ r . The second term in the integrand corresponds to cutting its (JV + /)-th line and 
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(N + l)-th column. By the same reason as before, the result of this integral is independent of 
the value of I, and it is equal to A^A^ -1 . 

Substituting the values of a^- and pij (eq.(4.1.2d)) in eq.(4.1.13), we finally get 



2N 

J 1=1 

_ oN-l 



for any value of /. 

Making the change of variables (4.1.3) in the integral H. 2 {1, k), it becomes: 



1 

H 2 (l,k)= QviOPou[ Qv 2 iPiv' 2 x 

"1,^=0 

2N 

2JV J2 Vi Dij v'j 

x / II d Vidf)'i e'' J=1 vl lN +i v' v [n+i Vu 2 N+k vl' 2 N+k- (4.1.15a) 

J i=i 

Once the i-th line and the i-th column need to be cut simultaneously if integral (4.1.15a) is 
to be non-vanishing, we have two situations to consider: 
i)l ^ k. 

In this case, the non-zero terms of eq. (4. 1.15a) are: 



(4.1.14) 



Q^OPOV! 1v 2 lPlv 2 x 

^1,^2=0 

2N 

2N E fix D " V'j 

x / II dVidVi e 1 -^ 1 ff VlN+l vl lN +i vl 2 N+k vl 2 N+k- (4.1.156) 

J i=i 

Since the result of the integrals in eq. (4. 1.15b) i independent of the space indices / and k 
(I 7^ k), only the total number of cuts within each sector v is relevant to their evaluation. For 
example, for v\ = and u 2 = 0, there are two cuts in the sector v = of matrix D. Thus the 
result of this integral is A^~ 2 A^. For v\ = and i/ 2 = 1, there is one cut in the sector v = 
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and another one in the sector v = 1 of matrix D; hence, the integral is equal to X 1 ~ A 2 ~ , 
the same result as for the case vi = 1 and z/ 2 = 0. 

Using the previous reasoning and the values of and pij given by eq.(4.1.2d), we get, 

H 2 (l,k) =2 N ~ 2 , for l,k=l,2,---,N and I ^ k. (4.1.15c) 

= k. 

In this case, we must take into account that rjj = fj 2 = 0. For / = k we have to compute 
all the possible permutations of the product v'iViV'n+iV'n+v The om y non-null terms in 
eq. (4. 1.15a) are: 





Po^i 


Qi/ 2 i 


Pl^2 




Vl 


Vi 


VN+l 


VN+l => 


QooPoo Qii Pi i 


Vl 


VN+l 


VN+l 


Vi =>- 


^ooPoi qnPio 


VN+l 


VN+l 


Vi 


Vi 


QioPoi Qoi 10 


VN+l 


Vi 


Vi 


VN+l 


qiopooqoipn 



(4.1.16) 



To write the multivariable Grassmann integrals as co-factors of matrix D, we have to 
rearrange the product to the pattern fj[r][ f]' N+l r}' N+l . Because of the permutations, the second 
and fourth term in (4.1.16) get a minus sign. 

By the same reason as explained in the calculation of integral the result of the 

integrals is independent of /, and it is equal to 

H 2 (l,l) =2 N ~\ for l = l,2,---,N. (4.1.17) 
In summary, we have that 



r 2N E Vi A]\ v., 

n 2 (l,k)= / Yl dVi^Vi e I,J=1 ViVi VN+k VN+k, 

J T 1 




iil = k (4.1.18) 
otherwise . 
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Finally, we calculate the integral Hs(l,k). Under the change of variables (4.1.3), it 
becomes: 



/ X 

2 



"2. "2 =° 



2N J2 fj'j Djj n'j 

x / Yl drjidfjj e J ' J=1 ^jv+i ^iv+fe-i vl 2 N+k v'^n+i+i- 

J i=i 

(4.1.19) 

The only non-null integrals in Tts(l, k) are the ones where I = k — 1. The terms on the 
r.h.s. of (4.1.19) that contribute to it are those where v\ = u[ and u 2 = v' 2 . 
Using the values of qij and pij given by eq.(4.1.2d), we get 

H 3 (l,l + l) =2 N ~ 2 , l = l,2,---,N. (4.1.20) 
Therefore, we have that, 



TC 3 (l,k)= / Yl dm&m e 7 " 7=1 ViVi+i VN+kVN+k-i- 

J 1=1 

27V " 2 ' + 1 (4.1.20a) 

From the results (4.1.14), (4.1.18) and (4.1.20a) we evaluate the terms < £(T, 0)£(T, 1) >, 
< £(<r,0) > CT ,< T-(t,0)T+(t,l) >,< £(T,0)W(1) > and < W(0)W(1) >. Substituting the 
results in (4.1.10), we finally get: 



(4.1.21) 

where Xb = —\9HbB. 
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The result (4.1.21) is valid for any value of iV space sites, and any set of values: 
(Eq, t, U, n) that characterize the fermionic system and its interaction with the external mag- 
netic field B. 



4.2. The Exact /3 3 Coefficient in d=(l+l) Hubbard Model 

In section 4.1 we used the simpler case /3 2 to exemplify, in great detail, how to use the 
results of Appendix A to calculate the coefficients of the expansion (2.4). In this section, we 
present the results of Tr[K 3 ] for any number iV of space sites. 

The expression of (3.12) for n = 3 is 



Tr[K^ 



» 67V 

=/n 

J 1=1 



drjidrji e 1 ' 



Vi A u VJ 



X 



x /C®(r/, 77; v = 0) K®% 77; v = 1) K®% 77; v = 2), 



(4.2.1a) 



where Ajj aze(the entries of the block-matrix A (eq.(2.13a). For n = 3, we have that 



xN 



\ 



(4.2.16) 



/ 



OWxtv Inxn -INxN 

\^NxN QvxAT IatxAT 

The expression of the Grassmann function JC®(fj, rj; v) for the case of the periodic unidi- 
mensional Hubbard model is given by eq.(2.1.3). 

The r.h.s. of eq.(4.2.1a) is expanded into a sum of 64 terms, each of which contains a 
Grassmann multivariable integral. Their evaluation will be attained by the same technique 
used in section 4.1 for the case n = 2. 

For n = 3, the matrices A ao ', a =f, j, are diagonalized by the matrix P (eq.(2.18a)), 



P = 



— InxN IjVxAf IjVxAf 

\ t N xN (l + e^^NxN (1 + e^)ljVxAT ) 
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(4.2.2a) 



the inverse of which is 



/ 



*N x N 



— 1 



NxN 



*-NxN 



\ 



(1 + e 3 jl^xiv l NxN {l + e3)t NxN 

\(l + e^)l NxN t NxN (1 + e^^N x N / 
For n = 3, the diagonalized matrix D (eq.(2.18)) is 



/21 NxN Qv 



D 



(4.2.26) 



(4.2.3) 



QvxTV \ 
QvxTV (l + e~)lAr X 7v QvxTV 
V QvxTV QvxJV (1 + e~^ ± )l NxN / 

The eigenvalues of the matrix A " ", a =f, |, are: 2, (l + e^) and (1 + e - s - ). Each eigenvalue 
is iV-fold degenerated. 

The change of variables (2.19), for the case n = 3 becomes, 



v'=0 



and 



VvN+l — ^2 rfv'N+l Qu'v 
i/'=0 



(4.2.4a) 



(4.2.46) 



with I = 1, 2, • • • , N, and 



1 (1 + eiP) (1 + e^ 1 ) 

Pv'u = I -1 1 1 

1 (1 + e—) (1 + e—) 



and 



(4.2.4c) 




2-TTi 

3 



-1 1 

g^v = ^ | (1 + e^) 1 (1 + e"' 

(l + e^) 1 (1 + eT) 

Using the symmetries studied in section 3, it is simple to show which terms on the r.h.s. 
of eq.(4.2.1a) are equal. Besides, taking the same steps as in section 4.1, we get that the terms 
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[13]: < £(0)£(l)T-(2) >, < 5(0)T-(l)T-(2) >, < £(0)T-(1)W(2) >, < £(0)W(l)T-(2) >, 

< r-(0)T-(l)T-(2) >, < r-(0)T-(l)T + (2) >, < r-(0)T-(l)W(2) > and 

< T-(0)W(1)W(2) > vanish identically. 

Taking into account the symmetries of section 3 and the fact that the mentioned terms 
are zero, we get that Tr[K 3 ] is equal to 



Tr[K 3 ] =< £(0)£(1)£(2) > +3 < £(0)£(1)W(2) > +6 < £(0)T-(l)T + (2) > + 
+ 3 < £(0)W(1)W(2) > +6 < r _ (0)T + (l)W(2) > + < W(0)W(1)W(2) > . 

(4.2.5) 

All the terms on the r.h.s. of eq. (4.2.5) can be written as the sum of multivariable 
Grassmann integrals. These integrals are reduced to five different types of integrals. We 
will not detail the calculation of these integrals, however, for it follows the steps presented in 
section 4.1 to derive the results of the integrals 7ii(l),Ti.2(l, k) and Tis{l, k). 

The five types of integrals that appear in Tr[K 3 ] and their respective results are [14]: 

1) 

g x {l)= Hdrndfj! e 1 '^ fhm = 2 N ~\ (4.2.6a) 

J T — 1 



1=1 



for / = 1,2,---,N. 
2) 



r 3^ Au r,j 

Q 2 (l,k)= / Y[ dr]idf}i e J - /=1 fji rji fj N+k r] N+k 



i=i 



3) 



k = l, 1 = 1,2, --^N 
2 N ~ 2 , k^l, l,k = l,2,---,N. 



n 3AT tj! A I3 rjj 

Ji(l,k)= I Y[ d Vidfji e 7 ' J=1 ViVi+i VN+k VN+k-i 

i=i 

2 N ~ 2 , k = l + l, 1 = 1,2, --^N 
0, k^l + 1. 
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(4.2.66) 



(4.2.6c) 



4) 



J2(h,k,h)= / Yl d Vidm e 7 " 7=1 Vh Vh Vn+i 2 Vn+i 2 +i V~2n+i 3 t]2n+i 3 -i 

J T 1 



1=1 



5) 



2 N ~ 3 , h ^ k, h + h and l s = l 2 + 1 

2 N ~ 2 , l 3 = h and h = h - 1 (4.2.6gQ 
0, for any other case. 



G3(h,h,k) = / II d Vidf]i e 7 - J=1 r/^ r] h fj N+ i 2 rj N+ i 2 fj 2N +i 3 V2N+i 3 

J i=i 

> 2 N ~* h = h + h, h = h± h, h = h± h (4.2.6e) 
2 , h = h = h- 



Using the results (4.2.6) to calculate the terms on the r.h.s. of eq. (4.2.5), after some 
algebraic manipulation, we get 

Tr[K 3 ] = 2 2N [N* {(E - ^ + \{E* - H?U + ^(E - ^)U 2 + ±U 3 ) + 
+ N 2 {^(E - ^) 3 + ^ (E - tfU +l(E - + 3(E - ^)t 2 + 

+ N(l(Eo - tfU + l(E - n)U 2 + l(E - v)t 2 - hj\ 2 B + \t 2 U + ^C/ 3 )] . 

(4.2.7) 



5. The Grand Canonical Partition Function for the One— Dimensional Hubbard 
Model Up to Order f3 3 

In eq.(2.4), we write the expansion of the grand canonical partition function in the high 
temperature limit in terms of Tr [K n ] , namely 
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Z(P^) = Tr(e^ K ) 

= Tr[I - (3K] + ^Tr[K 2 ] - ^Tr[K 3 ] + • • • , (5.1) 

In reference [9], we calculated Tr[l — K] for the unidimensional periodic Hubbard model 
for a chain with N space sites, where N can assume any arbitrary value, 

Tr[l - K] = 2 2N [1 - N(3((E - + j)] • (5.2) 
From expression (4.1.21), we have that 



Tr[K 2 ] = 2 2N [N 2 ((E - /i) 2 + i(Eo - y)U + *L) + 

_/V 3 
+ - ((Eo - n) 2 + (Eo -fi) + X 2 B + 2t 2 + ^U 2 )]. (5.3) 

In section 4.2, eq. (4.2.7), we got 



Tr[K 3 ] = 2 2 ^ [N 3 ((E - + \e - rfU + ^(E - ^)U 2 + ±U 3 ) + 

+ N 2 (l(Eo - fi) 3 + y (Eo - v) 2 U +\(Eo- fi)\l + S(E - ^)t 2 + 

+ l(Eo-,)U 2 + \t 2 U + \uX 2 B + l A U 3 ) + 

+ N(l{Eo - tfU + 3 -(Eo- »)U 2 + \(Eo- ^t 2 - \u\l + \t 2 U + |[/ 3 )] . 

(5.4) 

Substituting (5.2-4) in (5.1), we get the grand canonical partition function Z((3; n) of 
the unidimensional periodic Hubbard model. However, we will use the Helmholtz free energy 
W(/3; n) to calculate physical quantities. The relation between Z((3; /x) and W(/3; n) is given 
by 

W(/3;/i)=lnZ(/3; M ). (5.5) 

The Helmholtz free energy of the one-dimensional periodic Hubbard model, up to order 

(3\ is 
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2 In 2+ 



r/3 1 



+ (\(Eo ~ /i) 2 + \\% + \(E - n)U + y + ^U 2 )(3 2 - ((E - + |)/?] + 0(/? 4 ). 

(5.6) 

The coefficients of the terms (3, f3 2 and (3 s in the expansion of W(/3; /x) are exact. 

Once we have the Helmholtz free energy in the high temperature limit, we can derive 
any physical quantity in this limit. Here, we calculate only three quantities: 

i) average energy per site: < h >. 

The simplest way to derive the average energy per site from the Helmholtz free energy, 
is to scale the constants: (E ,t,U, Xb) — > (olEq, at, all, aXs) and substitute in eq.(2.1.2) to 
obtain W(/5; /x; a). 



<h>{(3) = 



1 5W(ft^a) 



/3iV 



da 



a = l 



(5.7) 



From eq.(5.6), we get that 



< h > ((3) = E + -+ 

+ [ ~ t 2 - ^U 2 - l -E U +\un- \e 2 + l -E 0l i - (3+ 
+ [ 3 -t 2 U - itV + ^E U 2 -\u 2 »- ±U\% + ±UE> - \uEop+ 

+ Yq U ^ + ^ + \ t2E ^ 2 + W 3 )- ( 5 - 8 ) 

ii) difference between average numbers of spin up and spin down particles per site: < n-f > 
- < ni >■ 

From the definition of the Helmholtz free energy (eq.(5.5)), we have that 
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< n T > (/3)- <rii> (/?) 
Up to order (3 2 , we get from eq.(5.6) that, 



(3N d\ B 



(5.9) 



A 



< n T > ((3) - < n; > (0) = '-f (U(3 + A) (3 + 0([3 3 ) . 



(5.10) 



Hi) average of the square of the magnetization per site: 



< ml > ((3) = -§ < (n lT -n u ) 2 > 



B 2 

,1 ,2 1 



dW(P;fi) +2 dW((3 ]f i) 



d/j, 



dU 



(5.11) 



where B is the external magnetic field. 
From eq.(5.6), we obtain that 



< m\ > ((3) = \g V B [I + \U[3+ 

+ [ ~ \uEo + \iiU - ±U 2 + ±\ 2 B - \e 2 + - ^ 2 }f3 2 ] + 0(f3 s ), 



(5.12) 



where g is the Lande's factor and \xb is the Bohr's magneton. 



6. Conclusions 

Using the grassmannian nature of the fermionic fields, we extend the results of our 
previous work [9] to higher order terms of the expansion of the grand canonical partition 
function for unidimensional self-interacting fermionic models in the high temperature limit. 

We are able to write the path integral for such models as a sum of co-factors of a suitable 
matrix, the entries of which are commuting quantities. This approach avoids ambiguities like 
the ones yielded by Hubbard-Stratonovich transformation [3] for fermionic path integral. 
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We apply the developed method to the Hubbard model in one-space dimension with 
periodic boundary condition. We get the exact coefficients of the (3— expansion (f3 = of 
the Helmholtz free energy up to order f3 3 for a ring with N space sites, where N is arbritrary. 

The expansion (2.4) of the operator e~^ K is valid for any value of (3. However, a sound 
physical meaning can only be granted to its lower-order terms if it is assured that the expansion 
converges. For finite values of the constants, there exists such a range: the high temperature 
limit. 

The calculation of the coefficients of (3 A and (3 5 terms of the grand canonical partition 
function for the one-dimensional Hubbard is to appear soon. 

The method presented here can be applied to any one-dimensional self-interacting fermi- 
onic model. We believe that this method can be extended to higher space dimensions. 
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Appendix A 

Moments of Gaussian Grassmann Mult i variable Integrals 



It is a known result [11] for a Grassmann algebra of dimension 2 2Ar , composed of the 
generators {771, • • • , rj N ; fj 1 , ■ ■ ■ , fj N }, that 

N 

I \drndfji e^ =1 =det(A), (A.l) 

' i=i 

where Aij are the entries of matrix A and are commuting quantities. 

We will show in this Appendix that the moments of integral (A.l) are co-factors of A. 

We first consider the case where we have one product f\\r\h in the integrand of the gaussian 
integral (A.l), that is, 

N 

/N ViAijijj 
Yldrjidfji fjirjk e^ =1 , (A2) 

i=l 

where /, k are fixed and 1 < I, k < N. 

Due to the fact that for all Grassmann generators we have: fjf = r/f = 0, i = 1, • • • , N, 
the only non-null terms in eq.(A.2) are the ones where the integrand has N products of the 
form: fji^j. Eq. (A. 2) becomes: 



M (l, k) = J JJ drjidra r}^ _ ^ A hjl ■ ■ ■ A lN _ ljN _ 1 x 

i=l V ii,---,ijv-i=l 

Ji,--',Jjv-i = 1 

X VilVjl Vi2 7 lj2 ' ' ' VlN-lVjN-l ) 

(A3) 

and the indices are such that i n 7^ I, n = 1, • • • , N — 1, and j n 7^ k,n= 1, • • • , N — 1. Once the 
product fji n r]j n is a commutative quantity, each term in the sum of (A. 3) appears (N — 1)1 
times. 

The (N — 1)1 distinct terms in (A. 3) can be generated by fixing one configuration for 
{ii, %2, • • ■ , ^at-i}, for example, we choose: \i\ = 1, ■ • ■ , ii-i = I — 1, i\ = I + 1, • • • , %n-i = N}, 
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and, taking all the terms coming from the sum over the indices j n ,n= 1, • • • , N—l. Therefore, 
M(l, k) becomes 



N N 



M ( l i k ) = / 11 dr li d Vi VlVk Al ^ ' ' ' A l-l,3i-i A l+l,3i ■ ' ■ A Nj N -i x 

i=l 3i,---Jn-i =1 



x mvh V2V32 ■ ■ ■ vi-iVji-iVi+iVji ■ ■ ■ VNVjN-i ■ i AA ) 

Renaming the variables: j t -> -> j z+2 , • • ■ Jn-i -> Jjv, we have that: 



r N N 

M(l,k)= Yl drudfji fjirik A iji • ' • A i-i,ji-i A i+i,ji+i- • ' A Njr* 

J i=i ji,"-,j"i_i=i 



J'l + 1,"-,JJV=1 

Defining the matrix B(Z, fc) as: 

and i, j = 1, 2, • • • , AT. 

Using the definition of matrix B(Z, A;), the expression of M(Z, fc) is re-written as: 



„ JV JV 

fc ) = / n s ■ ■ ■ B i-iji-i B iji ■ ■ ■ bn j^ x 

»■ 1 - 1 



x fkVji ■■■ViVh ■■■VNVjN- ( A -7) 
Integrating over fji, and using the definition of determinant [10, 15], we finally have that 



M{1, k) = detB 

= (-l) l+k A(l,k), (A.8) 

where A(l, k) is the minor determinant of matrix A, when the line / and the column k are 
deleted. M(Z, k) is the cofactor of matrix A. 
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Using an analogous procedure, we now consider the case of moments of the gaussian 
Grassmann multivariable integral when we have m products: 



ViiVk! VhVk 2 ■■■m m Vk m 

in the integrand of (A. 2), where m < N . 

Consider the fixed sets: L = {h,h, " • ,lm} an d K = {k\, k%, • • • , k m }. We define 
M(L, K) as 

N 

r N Yl ViAijVj 

M(L,K)= Y[d Vi dfj i fj ll rj kl ---fj lm r] krn e^ 1 , (A.9) 

J i=l 

and the products are ordered such that: h < fo < • • • < l m and k\ < ki < ■ ■ ■ < k m . 
Using an analogous reasoning, we obtain 



M(L, K) = J f[ d Vi dfji ' ' ' B i~hh-i B i,h 

*=1 Jir--,jjv = l 

x mVh ■ ■ ■ ViVji ■ ■ ■ VnVjn = detB(L, K), 

= ^l~j( l i+h-\ H m ) + (ki+k 2 -\ l- fc ™)y4^ X) 



B N j N x 



(A10) 



where the matrix B(L, K) is defined as: 
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B lJ {L,K) = { 



if % ^ h,- ■ ■ ,l n and j ^ h,---k r 
5 Ul 5 jkl , iii = h or j = k x 



(All) 



^ $il m Sjk m , if i = lm Or j = k m , 

and i,j = 1, 2, • • • , N. A(K, L) is the determinant of the matrix obtained from matrix A by 
deleting the lines: {Zi, 1%, • ■ • , / n } 5 and, the columns: {k\, ki, ■ • ■ , k n }. 

In summary, we can say that the effect of the presence of a product fjirjk within the 
integrand of the gaussian integral (A. 2), is to replace the line / of matrix A, Ay, by 5j k , and 
its column k, A ik , by Su. In their turn, the determinants of matrices B(L, if), eq.(A.lO), are 
easily written in terms of determinants of matrices of smaller dimension. Hence products of 
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Grassmann generators cut down the dimension of the matrices the determinant of which we 
are to calculate. 
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